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ABSTRACT 



We explore the amplification of magnetic seeds during the formation of the first stars and galaxies. During gravitational collapse, 
turbulence is created from accretion shocks, which may act to amplify weak magnetic fields in the protostellar cloud. Numerical 
simulations showed that such turbulence is sub-sonic in the first star-forming minihalos, and highly supersonic in the first galaxies 
with virial temperatures larger than 10 4 K. We investigate the magnetic field amplification during the collapse both for Kolmogorov 
and Burgers-type turbulence with a semi-analytic model that incorporates the effects of gravitational compression and small-scale 
dynamo amplification. We find that the magnetic field may be substantially amplified before the formation of a disk. On scales 
of 1/10 of the Jeans length, saturation occurs after ~ 10 s yr. Although the saturation behaviour of the small-scale dynamo is still 
somewhat uncertain, we expect a saturation field strength of the order ~ 10~ 7 n° 5 G in the first star-forming halos, with n the number 
density in cgs units. In the first galaxies with higher turbulent velocities, the magnetic field strength may be increased by an order 
of magnitude, and saturation may occur after 10 6 - 10 7 yr. In the Kolmogorov case, the magnetic field strength on the integral scale 
(i.e. the scale with most magnetic power) is higher due to the characteristic power-law indices, but the difference is less than a factor 
of 2 in the saturated phase. Our results thus indicate that the precise scaling of the turbulent velocity with length scale is of minor 
importance. They further imply that magnetic fields will be significantly enhanced before the formation of a protostellar disk, where 
they may change the fragmentation properties of the gas and the accretion rate. 

Key words, cosmology: dark ages, reionization, first stars - magnetic fields - Dynamo - Turbulence - Stars: Population III - formation 



1. Introduction 

The formation of the first stars is generally regarded as a 
well-defined problem, as the initial condition s at z ~ 100 
can b e derived accurately from CMB data (e .g. 
2009) using linear theory dBertschingett 



chem i stry is primordial and well-understood ( e.g. lAbeletall 
1997t iGalli & Pallal 119981: IStancil et al ' 
Yosh ida et al.L 120061; ISchleicher et al.L r2008bt iGlover & Savinl 
20091) . In addition, it is often assumed that magnetic fields are 



Komatsu et al. 
119981) and the 

1998; Omukai, 2001; 



not yet present and that the hydrodynamical equations are suffi- 
cient to describe the star formation process. 

This assumption is not neccessarily true. Indeed, a va- 
riety of mechanisms exist to create strong magnetic fields 
during inflatio n, the electroweak or the QCD phase transi- 
tion (see e.g. iGrasso & Rubinstein! 1200 ll for a review). By 
means of the inverse-cascade, the magnetic power of these 
fields may have been shifted to l arger scales in case of 
non-z e ro helici ty (Brandenb urg et all Il996t IChristensson et al.L 
l200lt iBanerjee & Jedamzikt 12004 . Strong primordial fields 
would have profound implications concerning the thermody- 
namics of the post-recombinatio n universe, reionizati o n and 
the formation of the first stars (ISethi & S ubramanianl I2005L 



Machi da et al .1120061 Tashiro & Sugiyamal 120061 ; iTashiro et all 
2006r. ISchleicher et all l2008ai l2009allbl) . " 

In this paper, we will however explore the limiting case 
in which extremely weak seed fields have been produced be- 
fore recombination. In such a case, the dominant contribution 
to the magnetic field strength comes from astrophysical pro- 
cesses after recombination. Cosmological MHD simulations in- 
cluding an approximate treatment of the Biermann battery term 
suggest that the Biermann battery could cre ate seed field o f 
the order 10 18 G in the IGM at z = 20 dXu et all |2008). 
Additiona l seed fields may be created v i a the Weibel instability 
in shocks (ISchlickeiser & Shuklal 120031 iMedvedev et al.ll2004l 
Lazar et all l2009h . The impor tance of dynamos in co smic sheets 
has early been recognized bv iPudritz & Silkl (Il989l) . The simu- 
lations of IXu et aL I d2Q08h run from cosmological scales to the 
protostellar collapse phase in a primordial minihalo. In such a 
situation, the following mechanisms are available to amplify the 
magnetic field: 

- gravitational compression of the magnetic field, 

- the small-scale turbulent dynamo which amplifies seed mag- 
netic fields already generated from cosmological processes, 

- large-scale dynamos in protostellar and galactic disks, 

- the magneto-rotational instability (MRI). 
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Gravitational compression under spherical symmetry leads to an 
increase of the magnetic field strength with n 2 ^ 3 , where n de- 
notes the number density of the gas. If the collapse proceeds 
preferentially along one axis, for instance because of rotation 
or strong magnetic fields, the scaling is closer to » 5 . In real- 
istic cases, often intermedia te values are found ( Ma chida et al.L 
2006; Banerie e et all 120081) . This ampl ification mechan ism has 
also been identified in the simulation of IXu et all d2008l) . 

Large-scale dynamos typically require the presence of a 
galac tic or protostellar disk and act on re latively long timescales 
(see lBrandenburg & Subrama nian. 2005, for a review). In such a 
disk, an exponential growth may also be obtained from t h e ma, 
netorotational instability (MRI, see iBalbus & HawlevL 1991 
which may seed other l arge-scale dyn a mos w ith turbulence 
dTan & Blackmanl l2004t ISilk & LangeH 120061) . However, the 
length scale of the fastest growing mode decreases for decreas- 
ing field strengths, and in the presence of a viscous cutoff length, 
such amplification may not be possible. As a result, a minimal 
fiel d strength is required to drive the MRI in a pr otostellar disk 
(see lTan & Blackmaril2T)c31lSilk & Langerll2006l for a detailed 
discussion). 

A critical condition for any dyn amo growth is that the idea l 
MHD approximation is applicable. iMaki & Susal (12004 120071) 
investigated this question using detailed models for magnetic en- 
ergy dissipation via Ohmic and ambipolar diffusion to show that 
the magnetic field is frozen into the gas unless it is very strong. 
An approximate fit to their results yields a critical field strength 
of about B < 10~ 5 (n/10 3 cm 4 ) ' 55 G. Due to the subtle effects of 
lithium chemistry, the ionization degree does not drop exponen- 
tially at densities of ~ 10 9 cirT 3 , but stays almo st constant with 
increa sing density. The more recent study by iGlover & Savinl 
(2009) finds even higher ionization degrees at these densities. 
This implies that the ideal MHD approximation can be used dur- 
ing the collapse phase to describe the interaction of magnetic 
fields with matter. 

Deviations from this behaviour may however occur on very 
small scales, where ambipolar and Ohmic diffusion become in- 
creasingly important. Estima t es based on the non -ideal MHD 
models of lPinto et ail d2008l) ; iPinto & Galhl d2008l) imply that, 
even if the magnetic field is in equipartition with the gas, am- 
bipolar diffusion is important only on scales 4 orders of magni- 
tudes smaller than the Jeans length, and Ohmic dissipation oc- 
curs only on even smaller scales. As Ohmic and ambipolar dif- 
fusion depend on the field strength itself, these scales will be 
significantly smaller for weaker magnetic fields, so that the ideal 
MHD approximation can be savely applied. This implies that the 
magnetic Reynolds number Rm = vl/rj varies strongly during 
the growth of the magnetic field, but always fulfills the condition 
Rmss> 1 and Pi~m = v/rj > 1 . Detailed calculations concerning 
the ambipolar and Ohmic diffusion scales will be presented in 
a companion paper, in which we make use of the ionization de- 
gree obtained from a numerical simulation to provide an updated 
calculation of these scales for different field strengths. 

Gravitational col lapse is generally accompan ied by the pres- 
ence of turbulence dKlessen & Henne bellg^ 20091). which ma y 
for instance be described bv the theory of [kolmogorovl dl94ll) . 
Numerical simulations show that primordial star formation dur- 
ing the collapse phase occurs in a self-similar fashion, where 
the density profile at a given time is always well-des cribed by a 
Bonner-Ebert sphere with a flat central density core dAbel et al.L 
l2002t iBromm & Lo"eR 120031: lYoshida et all |2008). Similar re- 
sults have been found for present-day star formation (e.g. 
IBaneriee et all 2006; Ba neriee & Pudritzl f2007). The gas falling 
on these central cores leads to weak shocks up to Mach 1, which 



drive turbulence in the central density core. This is reflected in 
the inhomogeneities in the central core and the sub-K e pleria n 
angular momentum pr o files, as reported bv l Abel et all ([2002); 
Brom m & Loebl (l2003h : lYoshida et all |2008). Under such con- 
ditions, a strong tangled magnetic field may be generated already 
during the collapse pha se by the small -scale dynamo that was 
originally proposed by Kaza ntsevl (1 19681) . 

This dynamo provides a very generic means of amplifying 
magnetic fields and was also prop osed to be impor tant in the 
large-scale structure of the universe (Ryu et al., 2008). The field 
amplification is due to the random stretching and folding of the 
magnetic field lines in a turbulent random flow. In the kine- 
matic regime, the field grows typically on the eddy turnover 
time, f e d = l/v where I is a typical turbulent length scale and 
v is the turbulent velocity. In the context of galaxy formation, 
iBeck et all d 19941) proposed that it is the small-scale dynamo that 
produc es the seeds for galacti c large-scale dynamos. As pointed 



out by Arshakian e t alJ d2009l) . the small-scale dynamo can ef- 
fectively amplify weak seed magnetic fields by ~ 13 orders of 
magnitude on a timescale ~ 300 milli on years in the first galax- 
ies. Similar results were obtained bv Ide Souza & Opheri ((2010) 
from a direct solution of Kazantsev's equation. 

Capturing such dynamos in numerical simulations of proto- 
stellar collapse is extremely challenging, as it requires that the 
turbulent cascade is well-resolved and well-separated from the 
scale where MHD turbulence is numerically dissipated. State of 
the art numerical simulations of turbulence thus require a spa- 
tial re solution of at least 512 3 for a marginally resolved inertial 
range (Federra th et alll2008ll2010l) . Numerical simulations fol- 
lowing protostellar collapse, on the other hand, typically resolve 
the Jeans length and thus the high density region with about 16 
cells, rendering them unable to capture the potential amplifica- 
tion via the turbulent dynamo. Feder rath et alJ (120101) showed 
that at least 30 grid cells are require d to reso l ve turb ulent vor- 
tices. Simulations as performed by IXu et alJ ((2008) therefore 
cannot resolve the turbulence in the central core. 

In this paper, we explore the implications of the small- 
scale dynamo during the gravitational collapse phase within a 
semi-analytic framework, applied to the formation of the first 
stars and galaxies. We first review the theoretical background 
and numerical evidence for the small-scale dynamo in § 12 and 
present a set of analytic estimates. In § [3] we develop a quan- 
titative model concerning the small-scale dynamo action dur- 
ing the collapse process. This model is applied in § |4] both 
to minihalos and atomic cooling halos, taking into account the 
amount of turbulence that was found in numerical simulations. 
Phenomenological consequences from the generation of such 
magnetic fields are discussed in § [5] In a companion paper 
dSur et all l2010t) . we present numerical simulations confirm- 
ing the importance of dynamo amplification during gravitational 
collapse. 



2. The small-scale dynamo 

In this section, we introduce the small-scale dynamo by sketch- 
ing its analytic derivation as well as numerical simulations that 
examined its efficiency, to make these results accessible to a 
broader community. This will be combined with a summary on 
the most important results and a first estimates concerning its 
importance in the first galaxies in § 12.31 Readers that are already 
familiar with dynamo theory or only interested in the most im- 
portant results may thus directly proceed from § 12.31 
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2.1. Analytical arguments 

The small-scale dynamo was introduced by Kaz antsevl 
(fl968h in an analytic framework, w hich was improved 
subse q uently in various works (e . g. [S ubramanian, 1998, 
1999; IBrandenburg & Subramanianl l2005l lArshakian et ail 
2009). For a detailed review we particularly recommend 



Brandenburg & Subramanianl d2005l) . as we can sketch here only 



the main steps in deriving the main properties of this dynamo. 
The induction equation of the magnetic field B is given as 



dB 

— = Vx(vxB-i;VxB) 
ot 



(1) 



where v is the velocity field of the fluid and rj the Ohmic or 
ambipolar resistivity. The velocity field is typically decomposed 
into a stochastic field \j and a drift component Vd- The latter de- 
scribes the change of the fluid velocity due to the Lorentz force. 
Such a decomposition is reasonable as long as the kinetic energy 
dominates over the magnetic energy. As we will see, this is al- 
ways the case in the regime where the small-scale dynamo is op- 
erational. In principle, the drift velocity depends on the complex 
history of field lines that a certain fluid element has seen during 
its history. To make the problem analytically tractable, however, 
it is generally assumed that most of the history averages out due 
to the tangledness of the magnetic field. Then, the drift velocity 
points at least approximately into the direction of th e instanta- 
neous Lorentz force . It is generally approximated as ( Kazantsevl 



1968; Subramanian, 1998; Brandenburg & Subramanian, 2005) 



vd = -;— [(V xB)xfi], 

4-np 



(2) 



where t is the typical response time of the magnetic field 
and p the gas density. The stochastic field component \j, on 
the other hand, is assumed to be an isotropic, homogeneous, 
Gaussian random velocity field with zero mean. We further 
adopt the Markovian ap proximation, assuming that its correla- 
tion function is given as ( Kazantsevl [l 9681: [Subramanian, 1998; 
iBrandenburg & SubramaniaiJ 20051) 



{v i T (x,t)v i T (y,s)) = T i K r )5(t- S ), 



(3) 



with r — \x - y\ and 6 denoting the delta distribution function. 
The matrix T'i can be decomposed in a longitudinal part Ti(r), 
a transverse part T^{r) and a helical part C(r). In the presence of 
homogeneous, isotropic and Gaussian turbulence, a weak seed 
magnetic field will be dragged with the fluid and obey the same 
properties. We denote the equal-time two-point correlation func- 
tion of the magnetic field as M' j (r, t) = (B'(x, t)B j (y, t) and de- 
compose it into a longitudinal component Miir, t), a transverse 
component M^(r, t) and a helical component H(r, t). Defining 
r' = x' - y' and introducing the Kronecker de lta 5 1 ' and the 
totally antisymmetric tensor egk, one can write d Subramanianl 
1 19981 IBrandenburg & Su bramanianl 120051) 



M' J = M„ 



S' J - 



+ M L — + He ijk r* 



(4) 



As the magnetic field is divergence free, one can then show that 



1 dr 2 M L 
2r dr 



(5) 



ISubramanianl 1998; Brandenburg & S ubra maiiianl 12005 ) : 



8M L 2 d I 4 8M L \ 
~df = l*A rKN ^r-) + GML 
dH 1 d I A d 



Aa N H. 



(6) 
(7) 



These equations include a diffusion term with coefficient a^v = 
rj + T L (Q) - T L (r) + 2aM L (0, f). The coefficient a N = 2C(0) - 
2C(r) - 4aH(Q, t) describes the strength of the a effect, and the 
term involving G = —4{{T^lr)' + (rT£)' /r 2 ) describes the rapid 
generation of magnetic fluctuations by velocity shear. The prime 
denotes here the derivative with respect to r. 

We will now focus on the case of nonhelical turbulence, 
with C(r) - and H(r, t) = 0, as it is more straightforward to 
be treated analytically. The formalism can however be applied 
to helical turbulence. In that case, the turbulent dynamo acts 
in the same way on small-scales, but creates additional large- 
scale cor relations that can ac t as seed fields for a large-scale 
dynamo (ISubramanianl 1 19991) . We are looking for eigenmode 
solutions of Eq. © and thus make the ansatz *P(r)exp(2rf) = 
r 2 -\[knMl. Insertion in Eq. [6] leads to a time independent 
equation which formall y resembles the Schrodinger equation 
of quantum mechanics (|Kazantsevl 1 19681 [Subramanian. 1998; 
Brandenburg & Subramanian, 2005): 



d 2x ¥ 

= -k n — + U (ryV, 



(8) 



Under the assumption that the velocity field is locally 
divergence-free, one can show that Uo(r) = T'^+(2T' L /r)+K'^/2- 
{k' n ) 2 I {Akn)+2kn I r 2 . Of course, the assumption of a divergence- 
free velocity field is only approximately true for compressible 
gas. However, a turbulent velocity field can always be decom- 
posed into a divergence-free and a divergent component, and the 
divergent component is just neglected to simplify the analytic 
treatment. As discussed in the next subsection, numerical simu- 
lations that generally include both components still find the same 
results concerning the magnetic field amplification timescale. 

As a boundary condition, we adopt *P — > for r — > 0, oo. 
The solution of Eq. © depends on the stability of flow of the 
magnetized fluid, which is described by the magnetic Reynolds 
number Rm = vl/rj, which describes the relative importance of 
the interaction of the magnetic field with the velocity field com- 
pared to Ohmic or ambipolar dissipation. As discussed in the 
introduction, high magnetic Reynolds numbers can be expected 
during primordial collapse. For the solutions of Eq. [8] a criti- 
cal magnetic Reynolds number Rm cr ~ 60 exists for which the 
solution corresponds to a bound state with F — ( Kazantsev, 
1 1 9681: ISubramanianl 1 1 998L fl999T) . For Rm > Rm cr , T > modes 
can be excited, leading to an exponential growth of the magnetic 
field on the eddy turnover timescale // v. Such a magnetic field 
is then curved on the turbulent length scale /. The thickness of 
the flux ropes is more uncertain an d may depend on the critical 
Reynolds num ber for the dynamo ( Subramanianl 1 19981) or the 
resistive scales dSchekochihin et all 120041) . 

To obtain the value at which the magnetic field saturates, 
the evolution equation (0 must be solved for a stationary state. 
The analysis shows that the magnetic energy then corresponds 
to Rm a ' of the kinetic energy. The maximum magnetic field 
strength that can b e obtained in this way is thus given as 
(SubramaniaHll998j) 



From the induction equation ([]]), one can t hen derive the fol- , 

lowing evolution equation for Mi and H dKazantsevl 1 19681 B max = J47rpv 2 Rm^ 



-1/2 



(9) 
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We note that the saturation field strength is still somewhat un- 
certain; our main intention here is to show that saturation can be 
reached. We checked that this conclusion does not depend on the 
precise fraction for the saturation level adopted here. 

2.2. Confirmation from numerical simulations 

As pointed out above, the complexity of modeling magnetic 
fields, gas dynamics and their mutual interplay forces one to 
make simplifying assumptions in an analytic treatment. Thus, 
numerical simulations are required to test the analytic results. 
Evidence of small-scale dynamo action has so far been con- 
firmed in numerical simulations of forced MH D turbulence 
(lHaugen et all l2004bt ISchekochihin et all l2004h both for high 
as well as low magnetic Prandtl numbers. Such simulations were 
able to follow the evolution of the magnetic field even in the non- 
linear regime when the Lorentz forces become s trong enough 
to sat urate the dynamo. In the simulations of lHaugen et alj 
(2004b), exponential growth of the magnetic field commences 
once the magnetic Reynolds number exceeds a certain critical 
value which for Pr m = 1 is Rm ~ 35. This critical value is 
somewhat larger by a factor 3 compared to the Rm„. obtained 
from the Kazantsev model, which assumes a delta-correlated ve- 
locity field. Simulations starting with an initially coherent field 
find that the magnetic power spectrum follows a k 3 ^ 2 slope be- 
fore saturation sets in. The growth rate of the field increases 
with (Rm / Rm„ ) 1 12 for higher magnetic Reynolds number. Apart 
from this, the salient feature of the small-scale dynamo is the 
presence of highly intermittent and structured fields. Th i s prop - 
erty has been discussed in de t ail bvlSchekochihin et alj (|2004); 
Brandenburg & Subramanian (2005). The magnetic fields ap- 
pear to be in folds whose length is comparable to the box size 
and which reverses direction at the resistive scale. The grow- 
ing magnetic field appears highly intermittent in the sense that it 
has strong positive(negative) values only in a few places in the 
simulatio n domain. These res ults were confirmed at higher res- 
olution (|Haugen et all l2004ah and for different Mach numbers 
dHaugen et all l2004d) . making a robust case for the efficiency 
of the small-scale dynamo. 

2.3. Order-of-magnitude estimates 

In summary, previous investigations found that the small-scale 
dynamo 

- amplifies magnetic fields on the eddy timescale, 

- leads to a magnetic fields curved on scales up to the turbulent 
length scale, 

- and a thickness of the flux ropes either given from the re- 
sistive scales or the critical Reynolds number for dynamo 
amplification. 

As the analytic procedures discussed in § 12.11 are based on 
the simplifying assumption that the gas is incompressible, they 
cannot be directly applied to a situation in which the gas is col- 
lapsing. Rather than this, we will focus now on the main results 
that were obtained with analytical and numerical simulations and 
develop a new framework where we can apply them during the 
star formation process. 

In order to assess the potential importance of the small-scale 
dynamo in the first galaxies, we start with an order-of-magnitude 
estimate concerning the rapid build-up of magnetic fields. The 
amount of turbulence in t he first galaxies has bee n studied with 
numerical simulations by iGreif et a"fl d2008h and IWise & Abell 



1(T 6 
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Fig. 1. Estimates on the field strength that can be reached within 
a free-fall time as a function of scale. Shown are both the ac- 
tual field strength to which the magnetic field is amplified by 
the small-scale dynamo, as well as the maximum field strength 
at which dynamo action saturates. This saturation field strength 
increases continuously as a power-law until the scale where tur- 
bulence is injected, whereas the actual field strength has a maxi- 
mum between 10 and 100 pc. We consider both Kolmogorov and 
Burgers turbulence. 

(120071) : IWise et all (120081) . On spatial scales of ~ 200 pc, they 
find typical turbulent velocities of ~ 20 km/s. To assess how 
much such turbulence can amplify magnetic fields in the first 
galaxies, we need to calculate the number of e-foldings within a 
free-fall time that are available for magnetic field amplification. 
The free-fall timescale is given as f// = 1/ ^[pG, with p the total 
mass density and G Newton's constant. We assume that the tur- 
bulent velocity scales as v(l) oc f, with B = 1 /3 for Kolmogorov 
turbulence and (3-1/2 for Burgers turbulence, on scales smaller 
than the injection scale. Kolmogorov turbulence describes a sit- 
uation where the gas is incompressible and is thus applicable in 
particular for sub-sonic turbulence, whereas Burgers turbulence 
describes turbulence in the presence of supersonic shocks, where 
the gas is quite strongly compressed. 

We note that Burgers turbulence should be considered as a 
highly idealized situation, because numerical simulations show 
that turbulence even in the supersonic regime always con- 
sists of ro tational and compressional components of compa rable 
strength (lHaugen et al.Ll2004dlFederrafh et aTll2008ll2010h. and 
dedic ated studies bylKritsuk et all d2007l) . ISchmidt et alj d2009l) 
and lFederrath et al](l2010l) typically find power-laws in between 
the Burgers and Kolmogorov case. As a side note, we mention 
that purely irrotational tur bulence may not be able to drive the 
small-scale dynamo at all (Mee & Brandenburg, 2006), whereas 
realistic turbulence should always have rotational and irrota- 
tional components as discussed above. 

The eddy timescale is given as f e d = 1 1 v. Considering that the 
small-scale dynamo should not be amplified above the saturation 
value given in Eq. (0, the magnetic field strength after a free-fall 
time is given as 

fi = min^ expJy-^=J,fi max J, (10) 

with Bq denoting the strenth of the seed field. To estimate the po- 
tential impact of the small-scale dynamo, we evaluate Eq. dTOb 
for a baryonic number density of 0.02 cm 3 , corresponding to 
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the mean density at virialization. We take into account the con- 
tribution of dark m atte r to the total mass d e nsity. According t o 
iGreif et all (120081) and lWise & Abell d2007l) : IWise et alj f2008). 
we assume that the turbulence is injected on scales of ~ 200 pc 
with a velocity of ~ 20 km/s. On smaller scales, we expect 
that the turbulent velocity follows some typical scaling law as 
in Kolmogorov or Burgers turbulence. We explore these cases 
for definiteness. We adopt a seed field of Bo ~ 10~ 20 G, some- 
what b elow the value e xpected from a Biermann battery mech- 
anism (IXuet all [2008). The results are shown in Fig.Q] In case 
of Kolmogorov turbulence, the magnetic field is negligible on 
scales of a few hundred parsec, but increases rapidly towards 
smaller scales. A maximum is reached on scales of ~ 80 pc, 
where the field strength has increased by about 13 orders of 
magnitude. The maximum occurs because the magnetic field has 
reached the saturation level on this scale. On all smaller scales, 
the eddy timescale is smaller, decreasing as Z/v oc I 2 / 3 in the 
Kolmogorov case. The magnetic field is thus saturated at smaller 
scales. As the turbulent velocity decreases with decreasing scale, 
also the saturation value for the magnetic field strength decreases 
correspondingly, thus explaining the maximum at ~ 80 pc. 

For Burgers-type turbulence, the situation is similar. As the 
turbulent velocity decreases more rapidly with scale, the mag- 
netic field reaches the saturation level only on somewhat smaller 
scales of ~ 30 pc. Towards even smaller scales, it decreases 
more rapidly, as the turbulent velocity field and thus the satura- 
tion field strength decreases more rapidly with decreasing length 
scale. 

Of course, this estimate is only approximate, as the magnetic 
field should also be subject to gravitational compression, which 
reduces the coherence length but amplifies the field strength. In 
addition, the density will increase during the gravitational col- 
lapse, thus increasing the field strength at which the magnetic 
field saturates. We also expect that the scale on which the tur- 
bulence is injected changes during the collapse. If the magnetic 
field saturates, additional effects like turbulent decay may come 
into play. To account for these additional physics, we propose a 
more detailed model in § [3] 



3. Magnetic fields during turbulent collapse 

In this section we derive a model for the evolution of a turbulent 
gas cloud with initially weak magnetic fields. For this purpose, 
a model of turbulence and the evolution of the turbulent spec- 
trum during the gravitational collapse is required. In this sec- 
tion, we first develop such a m odel on the basis of analytical ar- 
guments and numerical results dWise & Abell 120071: IGreif et all 
2008- IWise et alll2008h . Based on this model, the growth of the 
magnetic field in the dynamo phase can be calculated. We also 
discuss the subsequent evolution in the saturated phase, in which 
the competition between turbulent decay, gravitational compres- 
sion and dynamo amplification governs the evolution of the mag- 
netic field. 

3.1. Turbulent collapse 

During the galaxy formation process, turbulence is generated by 
the release of gravitational energy and the infall of accreted gas 
on the inner, self-gravitating core. Under such conditions, the in- 
jection scale of turbulence is usu ally comparable to the size o f 
the system under consideration (iKlessen & Hennebellel [2009). 
Here, we focus on the central density core found in numeri- 
cal simulations of primordial star formation dAbel et all [2002; 



Bromm & Loebl 12003b [Yoshida et al., 2008). Its extent is com- 
parable to the Jeans length, which can be derived from the criti- 
cal mass required to make the gas cloud gravitationally unstable: 




Here, c s is the sound speed. This central core is no longer dom- 
inated by radial motions, but rather supported by turbulence 
generated in accretion shocks. Due to the continuous infall of 
gas, the turbulence will not decay, but is constantly replenished. 
Radial profiles of the turbulent velocity in the first galaxies in- 
dicate that there are some random fluctuations, but the order- 
of-magnitud e of the turbulent velocity does not change during 
the co llapse dWise & Abell l2007l iGreifetall I20081 IWise etall 
2008). For turbulence driven by accretion, one expects that the 
turbulent velocity is comparable to the infall velocity, and for 
a roughly isothermal density profile, the free-fall velocity is in- 
dependent of radius. We will therefore assume that, while the 
injection scale changes during the collapse, the injected velocity 
stays the same. On scales smaller than the size of the cloud, we 
expect that the turbulent velocity scales as a v oc f . As above, we 
will explore both Kolmogorov and Burgers-type turbulence. The 
evolution of the me an density with time is pre scribed as in the 
one-zon e models of lGlover & Savinl d2009l) and lSchleicher et all 
(2009c), which follow the evolution of primordial chemistry dur- 
ing the collapse phase. 

3.2. Evolution in the dynamo phase 

As discussed in the introduction, there are two mechanisms that 
can amplify the magnetic field in the collapse phase: The small- 
scale dynamo, and gravitational compression. In this subsection, 
we propose a Lagrangian framework that allows to treat these ef- 
fects simultaneously. We consider an array of fluctuation length 
scales /, with a corresponding array of field strength B,, and fol- 
low both the evolution of the magnetic field strength B, and 
its corresponding length scale Z, over time. For the initial field 
strength, we adopt a conservative value of 10~ 20 G on all scales. 

The initial length scales are chosen with a constant loga- 
rithmic spacing, with the largest length scale corresponding to 
the Jeans length Aj, and the smallest length scale correspond- 
ing to 10~ 4 /ly. This range was adopted to make sure that the 
integral scale always lies well-within the range of scales fol- 
lowed in the code. We have estimated th e resistive scale follow- 
ing the framework of iPinto et all d2008l) : iPinto & GaliU J2008), 
finding that it lies well below this value. We follow the evolution 
of 1000 length scales with their corresponding magnetic field 
strength. For each length scale U, we calculate the turbulent ve- 
locity v,(f) as described in the previous subsection. As the small- 
scale dynamo amplifies each of these fields on the eddy-turnover 
timescale, we describe this process by solving the following set 
of ordinary differential equations (ODEs): 

dBj(t) = Biify 
dt li/vi 

In addition, gravitational compression amplifies the magnetic 
field strength and contracts the corresponding length scale. 
Under spherically symmetric conditions, the magnetic field 
strength scales as B oc p 2 ^ 3 , whereas a more modest increase 
would be appropriate if the magnetic field could significantly 
distort the geometry and lead to collapse along one preferred di- 
rection. In our model, we consider an initial seed field which 
is too weak to be dynamically important. Further, even after 
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sufficient amplification, the small-scale dynamo will produce a 
highly tangled magnetic field, which will also not change the 
geometry. We therefore increase the magnetic field strength ac- 
cording to this scaling law after each timestep. In addition, the 
corresponding length scale will be compressed. For collapse un- 
der spherically symmetric conditions, we expect it to scale as 
I cc p 1 / 3 . During the turbulent collapse, the initial seed field will 
thus be continuously amplified, while its length scale is com- 
pressed by gravity. 

As the magnetic field within the collapsing cloud is always 
distorted by the gravitational collapse, the largest coherence 
length that can be achieved is always smaller than the Jeans 
length by some factor fd- The precise value of fd is uncertain, 
as this problem has not been examined with numerical simula- 
tions yet. Throughout this paper, we adopt a fiducial value of 
fd = 0.1. A variation of this value by a factor of a few will 
however not change our conclusions significantly. Indeed, even 
somewhat larger coherence lengths may be possible, but this 
needs to be explored with numerical simulations. When we de- 
termine the integral scale in the subsequent applications, i.e. the 
scale on which the magnetic field strength is largest, we will thus 
explicitly ensure that it cannot be larger than fyA j . When the sat- 
uration field strength is reached, we do no longer solve Eq. d 1 2b . 
but switch to the treatment described in the next subsection. 



3.3. Evolution in the saturated phase 

Once the magnetic field B, on scale /,■ is larger than the saturation 
value given in Eq. (O, it is no longer amplified by the small-scale 
dynamo. It is however still amplified by gravitational compres- 
sion according to p 2 ^ 3 , whereas the saturation field strength only 
increases as p 1 ^ 2 . It can thus in principle increase above the satu- 
ration level. In this case, it is however subject to turbulent decay. 

In MHD simulations of decaying turbulence without self- 
gravity, one typic ally finds a power-law decay of the total 
magnetic energy (Smith & Mac Lowl[l998t lBiskamp & Miiller, 
1999). For our model, however, we need a prescription for the 
magnetic field on a give n scale rather than the to tal magnetic 
energy. As shown e.g. by Subramania n et al.l d2006l) . the power- 
law behaviour of the total energy results from an exponential 
decay of the magnetic field on a given scale, which happens on 
the eddy-timescale. This is because the integral scale increases 
during the decay, as the energy on smaller scales dissipates more 
quickly. Due to the growth of the integral scale, the total energy 
thus decreases as a power-law. 

In the presence of gravity, these effects have not been in- 
vestigated explicitly. However, it seems likely that the picture 
of purely decaying turbulence will not hold under these circum- 
stances. In particular, it is likely that the integral scale cannot 
become larger than a fraction of the Jeans length, which we 
parametrized above as fd- Once saturation is reached, the inte- 
gral scale within the central core thus continuously decreases 
over time. 

We therefore do not expect a power-law behaviour as in the 
case of decaying turbulence without gravity. For the kinetic tur- 
bulence, in fact we expect that typical turbulent velocities do not 
change significantly, although their scale will continuously de- 
crease during the collapse. This is because the turbulence can be 
continuously replenished by accretion shocks. A magnetic field 
stronger than the saturation value can however decay. On a given 
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Fig. 2. The evolution of the magnetic field strength on the inte- 
gral scale during protostellar collapse in a typical minihalo with 
subsonic turbulent velocities shown as a function of the average 
density reached in the central core at a given time. We show 
cases corresponding to Kolmogorov- and Burgers-type turbu- 
lence. 



scale, this should occur on the eddy-timescale, similar as for dy- 
namo amplification. We thus describe it as 



dBiif) Bi(t) 



dt 



kin ' 



(13) 



For a magnetic field strength above the saturation value, its evo- 
lution is thus governed by Eq. ( fT3b . while it is governed by 
Eq. ( TTZb for weaker fields. As a consequence, the magnetic field 
thus tends to stay close to the saturation value, which increases 



1/2 



as p 



4. Application to the first stars and galaxies 

The model developed above is now well-suited to study the im- 
plications of the small-scale dynamo in the first star-forming sys- 
tems in the early universe. In this section, we apply it to mini- 
halos that may harbor the first stars, and atomic cooling halos 
which are often considered to be the first galaxies in the uni- 
verse. 



4.1. Evolution in primordial minihalos 

The first stars are suggested to form in primord ial minihalos with 
10 5 - 10 6 M Q between z ~ 30 and z - 15 (lAbel etaU l2002l; 



Bromm & LarsonL l2004t iGloveri l2005t lYoshida et 31.0 2008) 



These halos have initial temperatures of a few 1000 K, and an 
initial ionisation degree of ~ 2 x 10~ 4 , corre sponding to the relic 
electr on fraction left after recombination ( Sea ger et all 1999, 
2000). The main cooling mechanism is thus line emission from 
molecular hydrogen, which forms primarily due to the H chan- 
nel. 

The turbulent properties in the first star-forming miniha- 
los have not yet been studied to a satisfactory degree and thus 
provide a relevant uncertainty in this analysis. The presence 
of turbulence is however evident from the inhomogeneities in 
the central cloud cores found in a number o f simulations, and 
the sub-keplerian angular momentum profiles dAbel et all 12 002; 
iBromm & LoeR I2003L lYoshida et all l2008h . iTan & Blac kman 
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Fig. 3. Magnetic field amplification in the early collapse phase. 
We show the same quantities as in Fig. [2] with particular focus 
on the initial phase. 
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Fig. 5. The magnetic field as a function of scale when an average 
density of 10 crrr 3 is reached in the central core, in a typical 
minihalo, shown for Kolmogorov- and Burgers-type turbulence. 




Fig. 4. The evolution of the integral scale of the magnetic field 
as a function of the average density reached in the central core 
at a given time during protostellar collapse in a typical minihalo 
with subsonic turbulent velocities. We show cases corresponding 
to Kolmogorov- and Burgers-type turbulence and compare with 
the Jeans scale. 



Fig. 6. Thermal and magnetic pressure as a function of the av- 
erage density reached in the central core at a given time during 
protostellar collapse in a typical minihalo with subsonic turbu- 
lent velocities, for Kolmogorov- and Burgers-type turbulence. In 
both cases, thermal pressure dominates over the magnetic one. 



(2004) mentions the presence of shock-velocities comparable to 
the sound speed. 

In our turbulent collapse model, we will thus assume that a 
turbulent velocity equal to the sound speed is injected on the 
Jeans scale, leading to sub-sonic turbulence on the integrale 
scale of the magnetic field. We consider both Kolmogorov and 
Burgers-type turbulence. 

With the model obtained in § [3] we follow the evolution of 
the magnetic field spectrum, and plot the evolution of the mag- 
netic field strength on the integral scale, as well as the integral 
scale as a function of density, in Figs. [2] and [4] In the initial 
dynamo phase, the magnetic field strength increases rapidly by 
several orders of magnitude, as the eddy-turnover time on sub- 
pc scales is much smaller than the free-fall timescale. We show 
the evolution during this dynamo-phase in more detail in Fig. [3] 
The spectrum of the magnetic field during this early phase, at 
a density of 10 cm 4 , is shown in Fig. [5] As one can see, the 
integral scale is initially small, but increases up to ~ 3 pc at a 



density of ~ 10" cm . At that point, the magnetic field is in the 
saturated phase, where it increases roughly with n , while the 
integral scale decreases as the Jeans scale due to gravitational 
compression. We find the same qualitative behaviour for the 
Kolmogorov- and the Burgers-type turbulence, indicating that 
the precise scaling of turbulent velocity with length scale is of 
minor importance. In both cases, we find a rapid build-up phase 
and a subsequent saturation phase. For Burgers-type turbulence, 
the build-up is just slightly delayed, and the field saturates on a 
slightly lower level, as the typical velocity on the largest possible 
integral scale, fa A j, is somewhat decreased. 

In Fig. [6] we show the evolution of magnetic pressure for 
both cases. We find that the thermal pressure clearly dominates 
over the magnetic pressure in both cases. Magnetic fields created 
by the small-scale dynamo may thus not change the protostellar 
collapse phase significantly. However, they are strong enough 
to change fragmentation behaviour and binary formation in the 
disk phase, as discussed in § Magnetic effects therefore need 
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Fig. 7. The evolution of the magnetic field strength on the in- 
tegral scale during protostellar collapse in an atomic cooling 
halo with turbulent velocities of ~ 20 km/s on the injection 
scale shown as a function of the average density reached in the 
central core at a given time. We show cases corresponding to 
Kolmogorov- and Burgers-type turbulence. 



Fig. 9. Thermal and magnetic pressure as a function of the av- 
erage density reached in the central core at a given time during 
protostellar collapse in aan atomic cooling halo with turbulent 
velocities of ~ 20 km/s on the injection scale, for Kolmogorov- 
and Burgers-type turbulence. 




10° 10 2 10 4 10 6 10 8 10 10 10 12 
Number density n [cm" 3 ] 



Fig. 8. The evolution of the integral scale of the magnetic field 
as a function of the average density reached in the central core 
at a given time during protostellar collapse in an atomic cooling 
halo with turbulent velocities of ~ 20 km/s on the injection scale. 
We show cases corresponding to Kolmogorov- and Burgers-type 
turbulence and compare with the Jeans scale. 

to be considered for a correct assessment of the primordial initial 
mass function (IMF). 

4.2. Evolution in atomic cooling halos 

Atomic cooling halos are defined as systems with virial temper- 
atures of at least 10 4 K and have masses of M > 5 x 10 7 [(1 + 
z)/10r 3/2 M . Graviational infall can thus shock-heat the gas 
to ~ 10 4 K, leading to an increased ionisation degree up to 
10~ 2 . Under such conditions, atomic hydrogen becomes an im- 
portant cooling agent. In addition, the presence of an additional 
cold accretion mo de has been found in numerical simulations 
dGreif et all [2008), which results form accretion in dense fila- 
ments with enhanced fractions of molecular hydrogen. The gas 
can thus cool efficiently while it is accreted, thus staying at tem- 



peratures of a few hundred K. Due to the large column densi- 
ties, Lyman a photon s may be efficiently trapped in these sys- 
tems (Spaans & Silk], 120061) . The presence of a photodissociat- 
ing background may suppress the abundance and increase 
the gas temperatures (Omuka"I l2001h . Even the combination of 
these effects is however not able to prevent the gas temper ature 
from decreasing during the collapse (ISchleicher et all [201 Oh . 

As sho wn by iGreifet all (120081) and IWise & Abell d2007l) : 
IWise et all d2008h . the gas in these systems is subject to super- 
sonic turbulence. They find supersonic turbulence with Mach 
numbers up to 10, corresponding to turbulent velocities of ~ 
20 km/s. We will adopt this as a fiducial value for the turbu- 
lent velocity on the injection scale and again consider the conse- 
quences of Kolmogorov and Burgers-type turbulence. 

As shown in Fig. [7] the magnetic field strength increases 
even more rapidly in this case, due to the supersonic turbulent 
velocities in atomic cooling halos. The magnetic field saturates 
at densities still lower than 10 cm" 3 and evolves then with n ' 5 , 
like the saturation scale. As shown in Fig. [8] the integral scale 
also increases more rapidly, yielding values of ~ 10 pc at a 
density of ~ 3 cm" 3 . Towards higher densities, its evolution is 
dictated by gravitational compression. It thus scales as the ther- 
mal Jeans length. As for minihalos, the difference between a 
Kolmogorov and a Burgers-type spectrum is of minor impor- 
tance, as in both cases the saturation value is rapidly reached. 

In Fig. [9] we compare magnetic and thermal pressure in 
atomic cooling halos. Due to the supersonic turbulent veloci- 
ties, the saturation field strength is highly increased, and thus 
also the magnetic pressure. The difference between the thermal 
and the magnetic field strength varies between 0.5 and one or- 
der of magnitude. In the presence of additional coolants due to 
metal-enrichments, this difference will decrease further. As dis- 
cussed below in more detail, the magnetic field will be highly 
inhomogeneous in a three-dimensional configuration and may 
thus be dynamically important locally. The magnetic field may 
thus change the evolution in the collapse phase and also after 
the formation of a disk. Potential consequences are discussed in 
more detail below. 
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5. Discussion and conclusions 



5.2. Formation of the first galaxies 



We demonstrated in this paper that magnetic fields are gener- 
ated rapidly by the small-scale dynamo both in minihalos which 
form the first stars, as well as in atomic cooling halos which may 
harbor the first galaxies. In this section, we summarize the main 
results and discuss the main consequences and open questions 
which cannot be resolved in this semi-analytic framework. 



5. 1. Formation of the first stars 

The formation of the first stars was often examined by hy- 
drodynamical simul ations that neglected potential effects from 
magnetic fields (e.g IXbel et all l2002t iBromm & LarsonL l2004t 
lYoshida et al.L I2006I) . Magnetic fields have however been con- 
sidered to be importan t in the protostellar d i sk in presence 
of an efficient dynamo dTan & BlackmanL [2004; Silk & Langer, 
2006). The impact of a uniformly imposed magnetic field (Bo) 
in the primo rdial collapse has been studied by Machid a et alJ 
(2006. [2008h using direct numerical MHD simulations. In these 
simulations, a range of Bo = 10~ 6 - 10~ 9 G was used for differ- 
ent values of the angular velocity of rotation. A protostellar jet 
within a radius of about 0.02 AU was found to be launched for 
an initial Bo > 10 G at n = 10 3 cm 4 . In this paper, we showed 
that the small-scale dynamo leads to a magnetic field strength 



much la rger than the critical v alue of 10 9 G (n/10 3 cm 3 ) de- 
rived by Mach ida et al.l (120061) for the formation of jets and out- 
flows. H owever, as these mag netic fields are more tangled than 
those of Mac hida et al.l d2006l) . their results cannot be directly 
applied to ours, and additional numerical studies concerning tan- 
gled magnetic fields are required. The average Alfv'en velocity 
va - B/ ^Jp is typically smaller than the sound speed, though 
it may dominate locally because of fluctuations in the magnetic 
field strength. 

In addition, our study clarifies that the magnetic field 
strength in the protostellar disk is much higher than previ- 
ously anticipated . The c ondi tions for the MRI, as formulated by 
Tan & Blackma3d2004|) and|Silk«fcLangei] (2006), are thus ful- 
filled. This leads to the presence of turbulence in the protostellar 
disk, which may drive a large-scale <mi dynamo in the presence 
of some kinetic helicity, making the field stronger and more co- 
heren t. But also the MRI itself may further amplify the magnetic 
field dBalbus & HawlevLfl99lh . 



This has important consequences for the fragmentation 
behaviour of the disk. Detailed numerical studies of the col- 
lapse of magnetised molecular cloud cores in the con t ext of 
present-day star format io n (e.g . IHennebelle & Tevssiei 
Hennebelle & Fromand, 120081: IHennebelle & Ciardil 



2008 



2009; 



Mell on & LiL [2009) indicate that even modest field strengths 
can suppress binary formation and strongly favour the formation 
of single stars. Jets and magnetic tower flows are very effective 
in transporting away angular momentum and thus change struc- 
ture and dynamics of the protostellar accretion disk. On the 
other hand, numerical simulations exploring the interaction of 
turbulence generated by the MRI with gravitational instabilities 
indicate the excitation of additional modes and an effective 
reduction of the accretion rate, as wel l as the broadening o f 
spiral arms by the MRI turbulence (Fromang et al., 2004). 
Dedicated numerical studies exploring the combination of such 
effects in a primordial accretion disk will thus be required to 
understand the full impact on the stellar masses. 



The formation of the first galaxies is currently subject to much 
larger uncertainties than the formation of the first stars. This is 
because the initial conditions are not completely clear and the 
amount of metal enrichtment is not fully understood. The pres- 
ence of superso nic turbulence has how ever been convinc ingly 
demonstrated bv lGreifetal] d2008l) and lWise etail (2008) with 
cosmological simulations encorporating hydrodynamics and pri- 
mordial chemistry. Additional physics like supernova feedback 
may just enhance the amount of turbulence found there. The 
small-scale dynamo is found to be extremely efficient under 
these conditions and may magnetise the material during the col- 
lapse, with the magnetic pressure only half an order of magni- 
tude below the thermal pressure on average. Locally, the mag- 
netic field may even dominate in some places, as it is expected 
that highly inhomogeneous fields are gen erated from the small- 
scal e dynamo. Indeed, as discussed by [Subramanian (1999) 
and [Brandenburg & Subramanian ( 2005J), the magnetic field is 
highly inhomogeneous, reaching equipartition in about 10% of 
the volume. The saturation field strength we adopted above re- 
sults from a spatial average over the different local values. We 
thus expe ct fluctuations of the field strength by at least a fac- 
tor of 10 dWang & Abell 120091: iDubois & Tevssierl l2009h . Due 
to the increase of the integral scale, the magnetic field becomes 
more coherent in these systems, making a stronger case for the 
putative presence of jets and outflows. 

As in the case of minihalos, the formation of a disk may lead 
to the presence of an ato dynamo that makes the magnetic field 
more coherent on disk scales. It may similarly play a role by 
making angular-momentum transport more efficient and thus re- 
ducing the amount of fragmentation and suppressing binary for- 
mation. As discussed above, the magnetic pressure may locally 
dominate over the thermal pressure. In this case, the magnetic 
Jeans mass sets the critical scale for fragmentation. For small- 
scale turbulent fields, it is defined in analogy to the thermal Jeans 
mass as 



M 



J,B 



2Mp 



0.2 km/s 



(l0 3 cm" 3 ) 



■1/2 



s 3 

2 ' 



P 



(14) 



Here, the Alfv'en speed va = Bj -\j4np replaces the sound speed 
c s , as magnetic pressure support propagates with the Alfv'en 
speed. With 

VA = 2 .o km/s (^(-l-p, (15) 

the Alfv'en speed in the saturation phase is thus larger than or 
comparable to the speed of sound. 

The presence of such fields thus provides additional 
stability during the formation of intermediate-mass black 
holes, wh ich are often considered to form in such sys- 
tems (e.g. lEisenstein & Loebl |1995|, iKoushiappas et al., 2004; 
Begelman ^all 120061: ISnaans & Silk, 20061: IShang et all 12009: 
ISchleicher et all 120101) . Thus, even if fragmentation cannot be 
totally avoided in hydrodynamical simulations, the presence of 
magnetic fields may still give rise to larger seed masses. The de- 
tailed consequences however need to be assessed with number- 
ical simulations. Additional open questions concern the further 
evolution of the magnetic field on large r scales and the build-u p 
of galactic-scale fields, as discussed by Arshakian et al. (2009). 
Their model for the magnetic -field evolution in galaxies yields 
a number of predictions which can be tested with future radio 
facilities such as the skaEL which can thus constrain the forma- 

1 http://www.skatelescope.org/ 



10 



Schleicher et al.: Small-scale dynamo action during the formation of the first stars and galaxies. I. The ideal MHD limit. 



tion mechanisms of the small- and large-scales magnetic fields in 
new born and young galaxies. In this respect, cosmological sim- 
ulations that include an approximate treatment for the mean-fiel d 
induction equation, as performed by Dubois & Tevssier (2009), 
will be very important. 

An additional issue that needs to be addressed is the role of 
magnetic helicity. Magnetic helicity is a conserved quantity and 
affects magnetic field generation and decay. In the presence of 
helical fields, th e small-scale dynamo also creates correlations 
on larger scales ( Subramanian, Il99l [T9991) . The decay law for 
helical fields was derived by lHatoril d 1984b . It is independent 
from the large scale part of the spectrum (i.e. scales above the 
integral scale) and is generally less efficient due to helicity con- 
servation. During such decay, magnetic power is shifted from 
small to large scales via a n inverse cascade, thus in creasing the 
typical coherence length dChristensson et all 1200 lb . Our paper 
was conservative in the sense that we assumed magnetic fields 
with zero helicity, yielding a lower limit on the integral scale. 
In the presence of helicity, magnetic fields may be coherent on 
larger scales, making it more straightforward to drive large-scale 
jets and outflows. To assess this issue, the turbulent properties of 
the first galaxies need to be analyzed and understood in further 
detail. 

5.3. Further discussion 

Based on the estimates performed for this paper, it seems likely 
that the small-scale dynamo will be very efficient during the for- 
mation of the first stars and galaxies. The epoch of first star for- 
mation may thus also be the epoch where the first strong mag- 
netic fields formed in the universe. This may be important for our 
understanding of primordial star formation. We further speculate 
that this mechanism may not only apply to the very first galax- 
ies, but that the formation of any gravitationally bound structures 
lead to a sufficient amount of accretion-driven turbulence to am- 
plify magnetic fields. We plan to investigate this proposition fur- 
ther with numerical simulations. 
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